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^ . Abstract 

0^ , We study the possibility that the approximate time shift symmetry during inflation is promoted to 
the full invariance under time reparametrization t — > t(t), or equivalently under field redefinition of 
the inflaton (f> — > 4>(4>). The symmetry allows only two operators at leading order in derivatives, so 
that all n-point functions of scalar perturbations are fixed in terms of the power spectrum normaliza- 
tion and the speed of sound. During inflation the decaying mode only decays as l/o and this opens 
up the possibility to violate some of the consistency relations in the squeezed limit, although this 
violation is suppressed by the (small) breaking of the field reparametrization symmetry. In particular 
one can get terms in the 3-point function that are only suppressed by in the squeezed limit 

£f~) , ki — > compared to the local shape. 
00 

O 

Introduction 

o 

CN ' The approximate scale-invariance of correlation functions produced by inflation is due to the dilation 
isometry of de Sitter space combined with the approximate symmetry of the inflaton dynamics under 
time translation [1] 

X ■ t->t = t + const . (1) 

& 



In this paper we want to explore the possibility that this symmetry is promoted to the full time 
reparametrization invariance 

t -> t(t) . (2) 

Of course this symmetry can be a good approximation only during inflation while it must be eventu- 
ally broken, similarly to what happens with the standard symmetry ([I]), at the end of inflation, when 
reheating takes place. This symmetry has recently been studied in the context of Hofava gravity and 
its healthy extensions [2 |3l S] ■ In these references the scalar mode describing the preferred foliation 
has been dubbed 'khronon'. See [5j El O [H] for other possible connections between Hofava gravity 
and the creation of primordial curvature perturbations. 
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We will see that, once this symmetry is enforced, the inflationary dynamics becomes very con- 
strained and unconventional. In particular three features are worth stressing. 

1. All correlation functions of £ are fixed, at the lowest order in derivatives, by only two coef- 
ficients, which can be written in terms of the normalization of the power spectrum and the 
speed of sound of perturbations. This is in contrast with the general case, where at any order 
in perturbations one can write new operators. 

2. During inflation the mode wavefunctions have the same form as in Minkowski. This apparently 
suggests the lack of a proper production of scalar perturbations. However, as we will argue 
below, this is not true if one considers the inevitable transition to a phase in which the time- 
reparametrization symmetry is broken. 

3. The above feature leaves an interesting signature in the correlation functions of the model. 
Indeed, the "decaying" mode decays much slower than in the conventional case (as 1/a instead 
of 1/a 3 ). This has remarkable consequences for the squeezed limits of correlation functions: 
the standard single-field theorems hold, but only at first order in the momentum of the long 
mode. One finds corrections at first order and, in particular, one has a l/k\ behaviour of 
the 3-point function in the squeezed limit. Unfortunately, these effects are very suppressed 
and totally unobservable. Indeed, the field redefinition symmetry itself is such that a time- 
dependent background wave, which would violate the consistency relations, can be removed 
and set to zero. Therefore, these effects are not there in the limit of exact field redefinition 
symmetry and they will only appear once we consider the small breaking of the symmetry. 

Section [2] describes the construction of the action compatible with the t — > t(t) symmetry. The power 
spectrum is studied in Section [3j with some details left to the two Appendices. The 3- and 4- point 
functions are discussed respectively in Section 2] and \5\ while conclusions are drawn in Section [6) 



2 Derivation of the action 

We want to write an inflaton action in which the usual (approximate) symmetry <p — > (j) + c is 
promoted to the full invariance under field redefinition eft — > (f>((p). We are going to assume an 
exact de Sitter metric and take the decoupling limit Mp — > oo, in which the dynamics of the scalar 
perturbations can be studied without considering the mixing with gravity. We will check the validity 
of this approximation in Appendix [Bj The time dependent inflaton background defines a foliation 
and in the presence of <f> reparameterization invariance, the only invariant object is the 4-vector 
perpendicular to the foliation [3] 

a "* (3) 



which is indeed invariant under (p — > (p(<j)). Notice that we are requiring a non-zero time-like d^cf). At 
low energy the operators with the smallest number of derivatives will dominate. It is straightforward 
to realize that it is not possible to write an operator with a single derivative. With two derivatives 
we have 

(V^) 2 ; V M u"V^; V^Vu,; u»u u V „u p V v u p . (4) 
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The first two are the same by integration by parts (this is true in the de Sitter limit where the 
Riemann tensor is proportional to the metric). Another constraint comes from the fact that u M is 
hypersurface-orthogonal, so that the Frobenius theorem implies 

V M u v V^ = W + u^u v V^u p V u u p . (5) 

We are thus left with two independent operators. The action to lowest order in derivative — and any 
order in u^ 1 — can thus be written as 

S = \J d A x^j (M 2 Pl R -2k- Ml (V^ - mf + Mlu^vrV^UpVvUp) , (6) 

where M a and M\ are the two parameters of our model, besides the vacuum energy A which is 
driving inflation. We subtracted 2>H from the term proportional to M 2 to reabsorb its contribution 
to the vacuum energy in A (notice that the cross term oc V M u^ is a total derivative). This action 
gives, at lowest order in derivatives, all the n-point functions and it will be the starting point for our 
calculations below. 

Another equivalent way to describe the model is by following the general construction of [9]. 
Any inflation model can be described in terms of the metric, in the gauge in which the infiaton 
perturbations are set to zero. One has to write operators invariant under time-dependent space 
diffeomorphisms and (approximately) invariant under time translations [9] 

%i — > xi (x, t) ; t — > t = t + const. (7) 

Here we promote the symmetry of the inflationary action to [21 H] 

Xi -> Xi (x, t) ; t t(t) ; (8) 

the symmetry (f> — > <j)((j)) becomes invariance under time reparametrization, as in this gauge constant 
time surfaces coincide with the ones at constant infiaton. Notice that the time reparametrization 
symmetry forbids to write operators with g 00 , which are otherwise allowed by the symmetries ((?]). 
The action ([6]) can be written geometrically as 

S = — £ d 3 x AtVhN [R^ + KijK ij - \{K - 3H) 2 + aa;aM , (9) 

in terms of the ADM variables 

ds 2 = -N 2 dt 2 + h lj (dx i + N i dt)(dx j + N j dt) , Kij = — (hij - ViNj - VjN^j , (10) 

and aj = N^diN. Indeed in this gauge one has (u u \7 u u^) 2 = a^a* and (V^u^) 2 = K 2 , so that the 
equivalence of the two actions follows from the Gauss-Codazzi relation (R^ = R^ + K^K^ - K 2 
up to total covariant derivatives), with the identification (A — l)Mp = M 2 and aMp = M 2 . Notice 
that in this language there are four invariant operators with two derivatives: R( 3 ', KijK^ , K 2 and 
aid 1 . One can get rid of one with the Gauss-Codazzi relation, up to a redefinition of the Planck 
mass. We still have an additional operator compared to the previous description. Indeed R^ does 
not play any role in the decoupling limit. Even more: as it is clear when one changes to spatially flat 
gauge, where R^ only depends on tensor modes, this operator does not affect scalar perturbations 
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even departing from the decoupling limit, or at non- linear order. This operator changes the speed of 
sound of gravitational waves as it affects their spatial kinetic term, but its effect is anyway negligible 
unless its coefficient is of the order Mp (0). 

The reader may be puzzled by the fact that the symmetry under field redefinition is incompatible 
with the fact that inflation must end once a certain point in field space is reached. But the situation 
is not different from the case of the usual shift symmetry, which will be strongly broken at reheating. 
Also here we only assume the field redefinition symmetry to be a good approximation while inflation 
occurs and perturbations are generated. Notice that a strong breaking of the symmetry in a region of 
field space where reheating takes place will not spoil the symmetry somewhere else, as renormalization 
is local in field space. In other words, the symmetry is valid only in a limited range of field space 
and it is badly broken if one considers field redefinitions which are large enough to move the point 
out of the symmetric region 



3 Power spectrum 

To calculate the power spectrum we expand the action ([6]) at second order. Using the field redefinition 
symmetry we can assume to perturb around <j)Q = t, i.e. <j) (x, t) = t + tt (x, t), in an unperturbed de 
Sitter space, which is a good approximation in the decoupling limit. Notice that the action does not 
contain any term linear in tt, which implies that the unperturbed Universe we are expanding around 
is indeed a good solution. In conformal time we get(f| 

* = / ^dr, (^(cV) 2 ~ %d\)^ . (11) 

This result is pretty unconventional. First of all, compared with the usual free-field action, each 
term has two additional spatial derivatives. This is not worrisome as additional spatial derivatives 
do not introduce extra pathological degrees of freedom. Second, the action does not contain any 
T) dependence so that the field is not sensitive to the expansion of the Universe and behaves as in 
Minkowski space (though with a speed of sound which is, in general, different from the speed of light). 
Actually these two peculiarities in some sense cancel each other to give a scale-invariant spectrum. 
Indeed, we expect the mode functions to be of the Minkowski form, but with an additional factor of 
1/k because of the presence of the additional spatial derivatives. It is easy to get the wavefunctions 

which give a scale-invariant spectrum for tt at late times rj — >■ 0. The curvature perturbation £ is 
given by £ = — Hit so that 



1 Notice also that one cannot induce sizeable graviton non-gaussianities cranking up the coefficient of this 
operator: indeed its coefficient cannot become parametrically large compared to Mj>, as this would imply a 
superluminal propagation of tensor modes. 

2 The tt exchange may induce spatial non-locality when coupled to other fields, as discussed in [3]. This is 
not relevant for us as we are not interested in coupling with other particles in calculating primordial correlation 
functions. Spatial non-locality may be relevant in discussing the horizon and flatness problem. 
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Notice that the scale invariance of the power spectrum (and of higher-order correlation functions) 
can be justified by symmetry arguments [T|, since we are in exact de Sitter and the action is shift 
symmetric. Of course, a small tilt is induced if the field redefinition symmetry is slightly broken. 

The result is encouraging, but the reader may be suspicious of this derivation. After all, how is 
it possible that perturbations are created if the field behaves as in Minkowski space? To understand 
what happens, let us follow the classical dynamics of a given Fourier mode. Although it is not 
sensitive to the Hubble friction, its wavelength is stretched and it eventually becomes much longer 
than the Hubble radius. In this regime the frequency of the mode, which keeps on oscillating as in 
Minkowski, becomes much slower than the rate of the expansion of the Universe. This means that, 
on a Hubble timescale, the time-dependence of the mode can be neglected and, similarly, its space- 
dependence becomes very small in a Hubble patch. We conclude that the solution we are describing 
is an attractor since the effect of perturbations becomes smaller and smaller as time evolves. 

This also sheds light on the quantum mechanical behaviour. Although each Fourier mode effec- 
tively remains in Minkowski, hindering a classical interpretation, the fact that its frequency becomes 
much smaller than the rate of expansion means that one is sensitive only to tt and not to tt. It is like 
probing in a laboratory a harmonic oscillator with an experiment which is very short compared to 
the period of oscillation: it will only be sensitive to the probability distribution of the position, but 
not to the momentum. The difference with the standard situation in inflation is quantitative, but 
not qualitative. Usually the time dependence of the mode decays, compared with the Hubble rate, 
as a~ 3 and it can safely be neglected. Here it decays as a -1 . 

The same logic also implies another important result: the conservation of £ on super-Hubble 
scales during the reheating stage and later. Independently of the details of reheating, we can assume 
that it will be insensitive to tt which is exponentially small compared to tt. This means that locally we 
are following the same unperturbed solution, with £ describing the relative difference in expansion 
between different points. In Appendix [A] we verify these intuitive arguments in an explicit toy 
example. We will see, in the following Sections, that this slow decay of the decaying mode leaves 
some signature in the higher-order correlation functions, which is a quite distinctive feature of this 
model. 

Due to the field redefinition symmetry one can choose the background solution to be <po = —r] 
and perturb now around this background (j) = —rj + x- It is straightforward to express at linear order 
these perturbations in terms of the perturbations around cosmic time as % = 7r/a, and write the 
second order action in terms of x from equation (jlip 

S2<X) = I d 3 xd V a* (^(cY) 2 - ^(d\f - Ml<H\dx?) . (14) 

This is compatible with the results of |10| . where it was noted that the effective mass is that of a 
conformally coupled field; this is consistent with the fact that the equations of motion for the field 
are like in Minkowski. Moreover, note that this action gives a power spectrum for x which is still 
scale invariant (since % and tt are related simply by a function of time) but with an amplitude that 
decreases exponentially during inflation. Different choices for the background solution seem to give 
different answers for the power spectrum in spite of the field redefinition symmetry. The issue is 
settled by the fact that what is more closely related to observations is the curvature perturbation 
conserved outside of the horizon £ which is equal to tt up to a constant factor as computed in 
Appendix [Bj 
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4 The 3-point function 



As we saw in the previous Section, the power spectrum for the fluctuations is scale invariant and 
indistinguishable from the predictions of more conventional inflationary scenarios. Let us now study 
the 3-point correlation function which carries additional information. It is conventional to define 



(C^C^) = (2Tr) 3 5(h + k 2 + k 3 )F c (h,k 2 , k 3 ) 



(15) 



where translational invariance implies that the 3-point function must be proportional to the Dirac 
delta, and rotational invariance implies that the function F^, called the bispectrum of £, is a function 
only of the magnitude of the momenta. As discussed in the previous Section, the dilation isometry 
of de Sitter, together with the time shift symmetry implies that the bispectrum is a homogeneous 
function of degree —6. 

The 3-point function of the field perturbation tt can be computed using the in-in formalism. It 
is given by (see Ref. |llj ) 



1 e J -° 



,H int (V)dT,' 3 



7r (r)*)Te *J-oo+ie 



(16) 



where |0) is the Bunch-Davies vacuum, T ( and T) indicates time ordering (and anti time ordering), ry* 
indicates the time at which inflation ends, "Hint is the interaction Hamiltonian and e is an infinitesimal 
positive constant. At leading order only the cubic part of the interaction Hamiltionian contributes, 
and one can show that Hint = — Ant- Therefore one can use the third order piece of the Lagrangian 
to compute the three-point function using 



n* 



dr] 



7r 3 (?7*), / d 3 xC int (t,x) 



The interaction Lagrangian can be computed by expanding the action, equation 
We get, after several integrations by partjf], 



(17) 

to third order. 



S* 



d xdr] ■ 
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Ml(2diir'diird 2 Tr + n' d l d j nd i d j n) + M^(n'd i ir"d i ^ - 8^' djirdidjir) 



(20) 



This cubic action coincides with eq. (5.10) of [5] in the Minkowski limit. In order to compute the 
3-point function for £ we use the relation £ = —Hn, additional non- linear terms in this relation either 



3 The reader might be worried about the appearance of an interacting term in the action that contains 
explicitly a second time derivative acting on tt that cannot be removed by a partial integration. This would 
not be a problem for us as we are treating these higher derivative terms as small corrections to the free action. 
However, it was noted in [3] that this term can actually be reabsorbed by performing a field redefinition of 
the form 

7T = 7f + 7f7f', (18) 



leading to the following action 
1 



d^x dr] ■ 



M 2 X 



(it d 2 ird 2 it' - ^{drrfd 2 ^ + M 2 (^it{dit') 2 - jTt(dff') 2 - d^'d^d^ 



(19) 



This action produces the same three-point function, eq. ([21]) , since the field redefinition (|T5|) vanishes outside 
of the horizon. One expects this to be true at every order in perturbations since in the unitary gauge, eq. 
©, the number of degrees of freedom is fixed 0]. 
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involve higher derivatives, which vanish outside of the horizon, or are suppressed by slow-roll factors, 
see 112] , We thus obtain the following expression for the bispectrum: 



F{(ki, fa, fa 



p 



fa - 
Tsiklh ■ k 2 + fcf&i ■ fa) 



-fa ■ fa 



fa ■ fa 



Ml k\ 



+ cyclic perms. , (21) 



where fa = fa + fa + fa and P^ = H / (2M a M\) is the £ power spectrum, eq. (|13p . All the 
contributions but the last cannot be large and give an /nl ~ 1- The contribution from the last term 
on the other hand is proportional to M^/M| = Actually it is easy to estimate the effect of 

each operator of the cubic action (|2Up comparing them with the quadratic action when modes freeze 
(dt ~ H, di ~ H/c s ). The only operator that can give a parametrically large 3-point function is the 
last in eq. (j2~T|) . 

We find an interesting feature of the model: it gives a single potentially large shape with an 
amplitude controlled by a single parameter, namely (? s . We plot the shape of this contribution in 
figure [TJ 



r\r\F 




Figure 1: We 'plot the shape of the part proportional to of the 3-point function, equation ([2~T|) . 
as a function of the ratios between momenta r 2 = fa /fa and r% = fa /fa, multiplied by r\r\. The 
shape is normalized such that its amplitude is one at the equilateral point r 2 = r$ = 1. 

In order to understand the phenomenological implications of this result, let us first introduce a 
quantitative way of comparing bispectra. One defines the scalar product between two shapes as |13| 



„ „ Pi{fa,fa,fa)F 2 (fa,fa,fa) 
1 ^.Se, WW ' 
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where the sum is over values of the momenta that form a closed triangle. One can then define the 
"cosine" of two shapes in the following way 

cos(Fi, F 2 ) = -, Fl ' F2 . (23) 

If the cosine between two shapes is close to one, one expects the data to be unable to distinguish 
between the two; conversely, if the cosine between two shapes is very small, constraints on the 
amplitude of one of the shapes do not constrain the amplitude of the othejf]. 

In CMB data analysis, a crucial numerical boost is gained when looking for shapes which are 
factorizable, i.e. which can be written as monomials of k\, k 2 , and k^. The standard procedure 
when comparing a theoretical 3-point function with constraints from CMB data is then to look 
for a factorizable shape which has a large cosine with the shape generated by a given inflationary 
model. Such shapes are often termed templates, which can be expressed as linear combinations of 
the so-called local, equilateral and orthogonal templates (see refs. |13[ I14|) . The cosines of the shape 
depicted in figure Q] with these three standard templates are 

cos(F c ,F local ) = 0.17, (24) 

C0S(F C , -F equ ii atcra i) = 0.93 , (25) 

cos(F f ,F orthogonal ) = 0.49. (26) 

It is therefore a good approximation to take the shape as equilateral. Its amplitude can be read from 
the expression (f2~Tj) above 

5 1 

f ^ L = lOS^f • (27) 
The limits on the equilateral shape obtained from WMAP 7 data given in ref. |15| can be used to put 
bounds on c s : c s > 0.013 at the 95% confidence level. Notice that (the potentially large contribution 
to) the 3-point function has a fixed positive sign in this modejf|. This is the opposite of what happens 
in more conventional models with reduced speed of sound (^-inflation), where the operator which 
reduces the speed of sound gives oc —\j(? s . However in those models one has another operator 
which contributes to the 3-point function and can flip the sign of /JS ; in our case we have no freedom. 
It is worth stressing that, although the shape given in eq. (|2ip has a large overlap with the equilateral 
one, the result has no free parameter and thus represents a potential smoking gun of the model. 

If one calculates the contribution to the 3-point function of the second and third operator in the 
action (|20p . each of them, when taking the squeezed limit k\ <C k 2 ,k%, diverges like 1/k 2 (while the 
leading term discussed above goes as 1/ k\ ) . This seems to contradict the results of references |16| [XT] 
where it is shown that, ignoring small deviations from scale invariance, in single field inflationary 
models the squeezed limit of the 3-point function diverges like l/k\ in that squeezed limit. This is 
due to the fact that the proof given there relies on Ckiv) vanishing at least like k 2 /H 2 outside of the 
horizon, while equation (fT2|) shows that in the model we are studying here C^ 7 /) vanishes only like 
k/H. However, this 1/k 2 divergence cancels between the two operators as it is evident in (|2ip " 



For CMB applications, this statement can be made more precise by defining a "two-dimensional" cosine, 
which takes into account the geometry and the effect of the linear transfer functions, to get closer to what it 
is actually observed |13j . 

5 We are using the WMAP sign convention for /nl- 

6 We are indebted to Austin Joyce for pointing out an error in the first version of the paper. 
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This is not an accident, but a consequence of the field redefinition symmetry. Indeed, a homogeneous 
time-dependent background mode, which would lead to the violation of the consistency relations, 
can be redefined away using the symmetry. Our theory is invariant under 



This means that a time-dependent background e(t) can be removed, provided we also redefine the 
field 7T as above. This redefinition, however, is irrelevant at late times as e — > 0, so we can conclude 
that the background wave has no effect and thus we do not violate any consistency relation. We 
checked explicitly that terms obtained from the cubic action with ir — > ir + e(t) cancel with terms in 
the quadratic action (|lip that appear after ir — > ir + e(t)ir. 

Notice, however, that the slow decay of the modes still opens up the possibility to violate the 
consistency relations: it is enough to consider terms which violate the original field redefinition 
symmetry. Of course, we expect these terms to be suppressed but to be there anyway, as indicated 
by the small observed deviation from a scale invariant spectrum. For example, we can add to the 
action ((9]) the cubic operator in unitary gauge 



which is not invariant under the field redefinition symmetry. This operator starts at cubic order, so 
that it does not modify the mode evolution, and it gives a cubic term ir'(d 2 ir) 2 which will violate the 
consistency relations. One might hope that such a non-standard behavior leaves an observational 
signature for example in the scale-dependence of the halo bias |18[ [TU} I20j . However, even under 
the optimistic assumption that the effect is suppressed by a single power of slow-roll compared to 
the leading l/c 2 . term, the analysis of references |21[ [22] (though performed for a different model) 
indicates that the observation of this effect with such a small amplitude seems unfeasible with planned 
surveys. 



Given that all the correlation functions to leading order in derivatives are completely fixed by two 
coefficients, it is of some interest to look at the 4-point function. In this section we compute the 
4-point function focusing only on the leading contribution proportional to c~ 4 . This part of the 
4-point function is important since observationally it gives the most relevant contribution in the case 
of small c s . 

In order to compute the 4-point function we need the interaction Hamiltonian to fourth order, 
for which it is no longer true that Hint = —£i n ,t- Let us start by expanding the action ([6]) to fourth 
order 



t + vr -)• F(t + tt) ~ t + vr + e(t + vr) + . . . = t + vr + e(t) + e(t)ir + -e(i)7r 2 + • • • ( 28 ) 




(29) 



5 The 4-point function 




a 



) ■ (30) 
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4 4) = -Ml [ d 3 xd V (^(dTrfidir) 2 + -{dirfd^drf + ^(^VAtf 
J V 8 2tt 

1 16 3 

+-zAird i ird j irdid j ir H ^(A7r) 2 (<97r) 2 + -^tt' Airdiirdrf H ^-ttVAttAtt 

or 2ar a 2 2or 

+ ^"A7r(d7r) 2 - — Tr'AvrfSvr) 2 + 4$7r3 1 -7r$7r / 7 V J . (31) 
or a ar / 

The second and third order pieces of the action are given by equations (jlip and (|20p respectively. 
Throughout this section we will only keep those terms that give the largest contributions to the 
4-point function in the small c 2 case. As before, it is easy to estimate the amplitude of the 4-point 
function by comparing each term in the quartic action (|3Up and (|3ip with the kinetic terms (|11|) once 
the modes freeze (dt ~ H and <9j ~ H/c s ). The amplitude of the largest piece of the 4-point function 
can thus be estimated to be proportional to cj 4 , generated by those terms in the non-linear action 
which are proportional to M 2 and containing the highest number of spatial derivatives. Thus, we 
will keep only the last terms in equations (|20p and (|30p . 

As stressed above, in order to obtain the correct expression for the 4-point function one must 
explicitly compute the Hamiltoniarf] 

H(P,tt) = P^' — C(P,tt), where the generalized momentum 
(keeping only the most relevant pieces in the small c 2 case) is given by 

P = ^j = -M 2 dV + ^d\dw? ■ (32) 



cV a 2a 

A straightforward computation of the terms in the fourth order interaction Hamiltonian which could 
potentiallygenerate a 4-point function proportional to c~ 4 using equations (|30p and (|32p shows that 
it vanished_|. 

In principle, two types of diagrams can contribute to the 4-point function: exchange diagrams 
and contact diagrams. However, since the fourth order interaction Hamiltonian vanishes, there is no 
contact diagram and the vacuum expectation value for the four-point equal-time correlation function 
in momentum space is given only by the exchange diagrams. In the in-in formalism the 4-point 
function can be then computed as 



3 J —CO 

2Re(r dr/f dif (0|C 5l C fil C fi ,C£ 1 (»7)«S(i/)«S(if)|0> > ) , (33) 

V J — oo J — oo I 



7 Notice that the canonical variables satisfying the commutation relations after quantization are the field 7r 
and the generalized momentum P, and the Hamiltonian is a function of these variables. Wherever we write 
7r' in the explicit expression for the Hamiltonian, it should be understood as shorthand for the appropriate 
expression in terms of 7r and P. 

8 It is important to note that in the full fourth order action, equation (|30p . there are terms containing two 
time derivatives acting on the field 7r" in such a way that they cannot be eliminated by an integration by 
parts. Similarly to what we did in the case of the 3-point function, we could have removed these terms from 
the action by a suitable field redefinition 

TT = 7f + 7f7f' + 8~ 2 di \ -Tt'diTX 



a 

which vanishes outside the horizon and does not change the result for the correlation functions. 
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The third order interaction Hamiltonian can be read from equation (|20[) . We are interested in 
the piece that can give a contribution to the 4-point function proportional to 1/cf which, after an 
integration by parts, we write as 

The time integrations appearing in equation (|33|) can be performed using 

° dr' c _^ (p+fcl+fc2)T . f° dr" ci %L (p+k3+k4)T » = r2 M* 1 1 1 



a(r') cl{t") M\ 2p3 (p + h + A; 2 ) 2 (p + h + k 4 ) 2 ' 

(35) 

and 

dr' ii£A( fcl+fc2 _ p ) T / /" r dr" ji£A( p+fe3+fc4 ) T » _ 



r 2 M a 1 1 1 , nP + k 3 + k 4 



a(r') J-oo a(r") 



The 4-point function can then be computed using equations (J33J) to ([36 



^ m\ 2 P 3 (p + fc 3 + h) 2 U? + 2 fc 3 ' • (36) 



x | {p e (h ■ k 2 )(h ■ k) - 2p 3 kf(p- k 2 )(k 3 ■ h)) 
+ {2p 3 kl(k 1 ■ k 2 )(p- fc 4 ) - Ak\kl{p ■ k 2 )(p- k 4 )) 



WTk\rw~ 1 k {kt+2{p+kl+k2)) 



1 

n 

+ 23 perms. , (37) 



4( P + fe3 + fe4) 2 + W (fct + 2(P + fcl + fc2)) 



where p = fci + &2 . 

Eq. ()37p is suppressed in the squeezed limit and does not contribute to the consistency relation 

mSkAA^ - 4° (f) p{q)p(k)2 ■ (38) 

This is easy to understand since it receives contributions only from exchange diagrams: when d 2 ir' 
corresponds to the external leg going to zero it will be trivially suppressed by g 3 , when diir corresponds 
to the external leg going to zero it will be contracted with both the other external leg which has 
some momentum k and the internal leg with momentum —k — q, which cancel at leading order in q. 

We checked also at the quartic level that the symmetry (f28|) holdall and this will prevent any 
violation of the consistency relations. Violations are possible, like in the cubic case, if one considers 
quartic terms which do not respect the field redefinition symmetry. 



J For this check it is crucial to keep also the quadratic action and vary it with the last term of eq. ([28 
this term cannot be neglected, because in going to conformal time it also gives an He' contribution. 
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6 Conclusions and outlook 



Given the simplicity of single-field inflation, it is certainly worthwhile exploring all the possible 
symmetries that can be imposed on its dynamics and their phenomenological consequences. Here we 
have studied the implications of imposing an approximate field redefinition symmetry <p — > on 
the inflaton. The predictions are very sharp since — after fixing the normalization of the spectrum — 
all correlation functions depend only on the speed of sound c s and are somewhat unusual, as a 
consequence of the slow decay of the decaying mode during inflation. 

What we have studied represents another de Sitter limit of inflation, as inflation can (but need 
not) take place with the metric being exactly de Sitter. This parallels the case of ghost inflation |23j . 
while another example has been studied in [9|. Like in the case of ghost inflation, the dynamics that 
may be responsible for modification of gravity in the late Universe, can be applied to inflation. This 
is not surprising, as models of modification of gravity often involve a scalar which defines a preferred 
foliation of space-time. And this is exactly what we need for inflation. 

It is useful to think about this model as another corner of the EFT of inflation [9|. Starting 
from a general situation, the limit H — > kills the unitary gauge operator g 00 , and therefore the 
standard spatial kinetic term of the inflaton. This is the limit of ghost inflation [23], when the 
spatial kinetic term is given by higher order spatial derivatives (K 2 and K^K^), while a standard 
time kinetic term 7r 2 comes from the unitary gauge operator (g 00 + l) 2 . The symmetry that we 
discussed forbids any operator of the form (g 00 + l) n , so that also the time kinetic term is now given 
by the higher derivative operator N~ 2 (diN) 2 . Of course these are only limiting cases: intermediate 
regimes in which various operators are relevant may have interesting features. We leave this to future 
investigations. 

It is important to stress a relevant drawback of our model, i.e. its spatial non-locality: the Green 
function of ir shows instantaneous propagation of the signal as discussed in [24] , Most likely, this 
implies that our EFT cannot be embedded in a standard Lorentz invariant UV completion 10 !. This 
is similar to what happens in models of fc-inflation with superluminal speed of sound c s > 1. 
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A Evolution after the field redefinition invariant phase 



In this Appendix we want to verify our intuitive arguments of Section [3] in an explicit (toy) example. 
Let us add to the quadratic action (|lip a standard 2-derivative kinetic terrrF^I 



S= / d 3 xdr] 



(39) 



We need ^ < 1 for the kinetic term discussed in this paper to dominate at Hubble crossing. In 
this case /3 represents a small breaking of the field redefinition symmetry and the contribution of the 
kinetic term we added will become relevant when a mode is sufficiently long compared to the Hubble 
radius. What we want to check is that, up to corrections suppressed by (3, tt remains constant during 
the out-of-Hubble evolution, until the mode becomes long enough to be dominated by the standard 
kinetic term. This will imply that the correlation functions calculated in the paper are actually the 
ones observed at late times. The equation of motion is given by 

d\" - (3H 2 — [a VI - (8\ - (3a 2 H 2 d 2 ir) = . (40) 

Out of the Hubble radius, i.e. (k/aH) 2 <C 1, there are three regimes of different evolution. For 
^2/3 ^ (k/aH) 2 <C 1, the terms proportional to j3 are irrelevant and everything goes as discussed 
in the paper. The first term which becomes relevant is the Hubble friction and it is easy to realize 
that this is the only term one has to consider in addition to the original Lagrangian in the window 
/3 < [k/aH) 2 < /3 2 / 3 . Finally, in the regime (k/aH) 2 <C /?, only the terms proportional to j3 are 
relevant and tt behaves as in standard inflation. It is simple to follow the evolution from one phase 
to the other in the long wavelength limit. First of all notice that tt = const is a good solution in any 
phase and in the transition regions for a mode which is well outside the Hubble radius, i.e. in the 
k — > limit. This can be seen explicitly in the equation and follows from the general conservation 
of £ on super horizon scales (which within our approximations implies the conservation of tt as 
£ = —Htt, with constant H). Moreover, the velocity becomes irrelevant, tt <C Hit, before the terms 
proportional to j3 start playing any role, and this implies that tt can be neglected when matching to 
the next phase. There is no mode mixing and tt remains constant all along. It is easy to check this 
behaviour numerically. 

The same reasoning works if we allow f3 to be time-dependent, i.e. dependent on the background 
value of 0o = t. This describes the fact that the field redefinition symmetry will be badly broken 
at the end of inflation and f3 will become large. It is straightforward to check that also in this case 
tt = const is a good solution so that, for modes well ouside the Hubble radius, i.e. tt <C Htt, the 
field remains constant while the symmetry gets broken. Notice that the logic is exactly the same one 
uses in the case of standard inflation to justify the conservation of £ through the unknown reheating 
phase. As in that case we expect the same arguments to be valid non-linearly in the amplitude of £, 
so that each n-point function remains the same when out of the horizon. 



11 For simplicity we assume that the speed of sound of the kinetic term we added is the same as the one of 
the original terms. 
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B Constraints and the validity of the decoupling limit 



In the main text we have calculated everything in terms of ir, focussing on the decoupling limit, 
(i.e. neglecting its effects on the metric) and then converting the results in terms of £. The logic 
behind it is that we expect the corrections coming from the effect of tt on the metric to be subleading 
in 1/Mp, and therefore negligible when M? <C M|> and M% <C MJ>. However the model we are 
describing is sufficiently unconventional to warrant a check of this intuitive explanation. Let us 
calculate the power spectrum of £ directly in the £-gauge, i. e . setting to zero the tt perturbations. 

Starting from the action ([9]), we go through the standard procedure of solving the constraint 
equations and plug the solution back into the action. We use the ADM splitting of the metric (|10p . 
Defining N = 1 + 5N and N l = + diip, with diNj, = 0, the linearized constraint equations 
obtained by varying with respect to iV and N l respectively are given by 

+ H) (# ,_3 ( C- OT )) +# (^) + W^0 



di 



3M?\ M? , 

. (41) 



(5NH-C)[l + - T A +^d^ 



2MpJ 2M p 



We can now solve these equations at first order in M^/Mp and M?/Mp and plug the solutions back 
into the action. After some work, we obtain 

S C = /d 3 xd,(^(90 2 -^(9 2 C) 2 ) , (42) 

which is the action given in (fTT|) with tt = —Q/H as expected. The action above will contain 
additional terms suppressed by powers of M^/Mp and M^/Mp. 
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